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I. INTRODUCTION 

The fluctuations in Cosmic Microwave Background 
(CMB) contain an amazing amount of information about 
our universe. Detailed measurements of anisotropy in the 
CMB reveal global properties, constituents and history of 
the universe. In standard cosmology, CMB anisotropy is 
assumed to be statistically isotropic and Gaussian. Gaus- 
sianity implies that the statistical properties of the tem- 
perature field can be completely characterized in terms 
of its mean < AT >= 0, and auto-correlation function 
C(hi,h 2 ) =< AT(hi)AT(n 2 ) >, where ft = (6,(f>), is a 
unit vector on the sphere. The angular brackets < .. > 
denote ensemble expectation values, i.e, averages above 
are for all possible realizations of the field over a sphere. 
Since we have one CMB sky, that is just one out of 
all possible realizations, the ensemble expectation value 
C(hi, n 2 ) can be estimated in terms of sky averages only 
to a limited extent, depending on underlying symmetries 
in C{hi,fi2). Under the usual assumption of Statistical 
Isotropy (SI), implying essentially Einstein's cosmologi- 
cal principle for cosmological perturbations, the correla- 
tion function is invariant under rotations. It implies the 
correlation function C(fii,n2) — C(fi\.ri2) = C{&), can 
be readily estimated by averaging over all pairs of sky 
directions separated by an angle 9. 

Spherical harmonics form a basis of the vector space 
of complex functions on a sphere, making them a natural 
choice for expanding the temperature anisotropy field, 

AT(ft) =J2 a imYi m (n)- (1) 

lm 

Here AT is the temperature fluctuation around some av- 
erage temperature T. The complex quantities ai m are 
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drawn from a Gaussian distribution, related to the Gaus- 
sian temperature anisotropy as 

ai m = / dn A Y lm (h)AT(h) . (2) 

The condition for SI now takes the form of a diagonal 
covariance matrix, 

< a lim 1 &i 2m2 >— C[ 1 Sl 1 l 2 S mim2 . (3) 

Here C/ is the well known angular power spectrum. In 
the SI case, the angular power spectrum carries complete 
information about the Gaussian field, and the statistical 
expectation values of the temperature fluctuations are 
preserved under rotations in the sky. This property of 
CMB has been under scrutiny since the release of the first 
year of WMAP data. Tantalizing evidence for statistical 
isotropy violation in the WMAP data using a variety of 
statistical measures has also been claimed in recent lit- 
erature [IHij]. However, the origin of these 'deviations' 
from SI remains to be modeled adequately. These devia- 
tions could be either genuinely cosmological, or statistical 
coincidence, or residual foreground contamination, or, a 
systematic error in the experiment and the data process- 
ing. Hence, it is important to carry out a systematic 
study of SI violations using statistical measures within a 
unified, mathematically complete, framework. Moreover, 
it is important to develop several independent statistical 
measures to study SI violations that can capture differ- 
ent aspects of any measured violation and provide hints 
toward its origin. 

While testing a fundamental assumption, such as SI, 
is in itself a justifiable end, there are also strong theo- 
retical motivation to hunt for SI violations in CMB on 
large scales. Topologically compact spaces (34To| and 
anisotropic cosmological models QjHl3[, Cosmological 
magnetic fields generated during an early epoch of in- 
flation 14| can also lead to violation of SI [l5[ , are but a 
few examples. 

This paper focuses on linking measures of SI viola- 
tion to the reduced symmetries of the underlying cor- 
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relation patterns 1 in the CMB map or the correlation 
function. While we present illustrative examples of the 
symmetries from various mechanisms of SI violation, this 
paper does not concern itself with a study of specific 
mechanisms. We define, within the framework of Bipo- 
lar harmonic representation of CMB sky maps, a num- 
ber of observables that can be used to quantitatively test 
SI. We present a study of the properties of bipolar mea- 
sures as one systematically reduces the rotational symme- 
tries of the CMB correlations, as is expected in different 
theoretical scenarios. We recapitulate the bipolar har- 
monic representation and the definitions of a set of mea- 
surable quantities representing SI violation in section [TT1 
In Section Hill these observable measures are computed 
for different levels of residual rotational symmetries of 
CMB correlations. This provides a clear understanding 
of the underlying symmetries revealed through the differ- 
ent bipolar measures. Section llVI deals with bipolar for- 
malism measurables using Bianchi template. Section [V] 
summarizing conclusions and discussions is followed by 
appendices where details of the calculations leading to 
results are presented. 

II. BIPOLAR FORMALISM AND THE 
OBSERVABLE MEASURES 

Any deviations from SI introduces off-diagonal terms 
in the covariance matrix Eq. ([3]), thereby making Ci an 
inadequate quantity to characterize the statistical prop- 
erties of the temperature field [lj|- Under such a sit- 
uation Bipolar spherical harmonic expansion, proposed 
by Hajian and Souradeep [l7], Il9l - [2"3l |. proves to be the 
most general representation of the two point correlation 
function, where the angular power spectrum Ci is a sub- 
set of Bipolar spherical harmonic coefficients (BipoSH). 
Two point correlation function of CMB anisotropies can 
be expanded as 

C(n h n 2 ) = Yl A ifA Y h(ni)®Yh(n 2 )}LM , (4) 

h,l 2 ,L,M 

here Aj^f 1 are Bipolar Spherical Harmonic coefficients 
(BipoSH) 2 , |/i - l 2 \ < L < (h + l 2 ), mi + m 2 = M, 
and {Yi 1 (hi)®Yi 2 (77.2)} lm are Bipolar spherical harmon- 
ics [24]]. Bipolar spherical harmonics form an orthonor- 
mal basis onS 2 xS 2 , with transformation properties un- 
der rotations similar to spherical harmonics. The tensor 
product in harmonic space can be explicitly written using 
Clebsch-Gordan coefficients C^^j as, 

{Y h {hi) ®Y h {n 2 )} LM = 

E C hm 1 l 2 m 2 Y hm x (ni) Yi 2m2 (h 2 ) . (5) 
m\m 2 



In this paper, we use the term 'correlation patterns' to inter- 
changeably refer to SI violation 



A. Bipolar spherical harmonic coefficients - 
BipoSH 

BipoSH can be extracted by inverse transformation 
of Eq. dU, i.e., multiplying both sides of this equation 
by {Yi'^hx) ® Y[' 2 (fi2)}* L , M , , and using orthonormality of 
Bipolar spherical harmonics. Hence, given a real space 
correlation pattern BipoSH coefficients can be found us- 
ing 

Af£ = J dSh 1 dO th C{h 1 ,fh){Yi 1 (h 1 )9Y ta (rh)}l M . 

Since C{n\n 2 ) is symmetric under the exchange of h\ 
and h 2 , this gives rise to following symmetry properties 
of BipoSH: 

a LM _ ( 1 \li+l2 — L aLM 

Af t M = Afi M 6 L , 2k . fc = 0,l,2,3,....(6) 

Hence, Afj M exists for even L and vanishes otherwise. It 
was shown in [l?} that the Bipolar Spherical Harmonic 
(BipoSH) coefficients Aj^f^ are a linear combination of 
elements of the harmonic space covariance matrix includ- 
ing the off-diagonal elements that encode SI violation, 

^hh = E < a hmx a * 2 m 2 > (~ l) r ™ 2 Q^mi2 2 -™2- CO 
mi m 2 

When SI holds the covariance matrix is diagonal, Eq. ^ 
and Clebsch property (IH6|) . therefore 

Affi = (~l) h C h (2h + l^^Ao^/o, (8) 

implying that A^ contains all the information on the 
diagonal harmonic space covariance matrix given by C;. 

The well knownpower spectrum Ci thus forms a sub- 
space of BipoSH [Ijj]. Under SI, the only non-zero Bipo- 
lar spherical harmonic coefficient will be Af® (equivalent 
of Ci), all the rest must vanish. The violation of SI thus 
implies A®® are not sufficient to describe the field. Hence, 
BipoSH proves to be a better tool to test SI, as non-zero 

A,j, other then L — and M = 0, terms should confirm 

1 2 . . 
its violation. 

It is impossible to measure all Aff from just one CMB 
map because of cosmic variance. Thus we need to com- 
bine them in different ways to diagnose different aspects 
of SI violations. 



B. Bipolar power spectrum- BiPS 

The Bipolar Power Spectrum (BiPS) is a rotationally 
invariant, quadratic measure that can be constructed out 
of BipoSH coefficients [l?} • BiPS involves averaging over 
BipoSH that reduces cosmic variance in comparison to a 
single CMB map, however this does not erase all the SI 
signatures. BiPS is defined as 

h,h,M 
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For statistically isotropic models Kl = kqSlq, i.e., kl — 
V L > 0. Thus a breakdown of SI will imply non- 
zero components of BiPS. In real space, kl can be ex- 
pressed as 



( 2L + 1 



dfl &1 



J dR X L {R)C{Rh 1 ,Rh 2 ) 



(10) 



where C(Rh\, Rh 2 ) is the correlation function after ro- 
tating the coordinate system through an angle w (0 < 
uj < 7r), about an axis n(8,$). Rh\ and Rfi2 are the 
coordinates of the pixels h\ and h 2 in the rotated coordi- 
nate system. The rotation axis n, is characterized by two 
parameters (0 < 9 < n), and $ (0 < $ < 2tt). x L , is 
the trace of finite rotation matrix in LM-representation 
called the characteristic function, and it is invariant un- 
der rotation of coordinate system, 



x l {R) = Y j d mm{R)- 



(11) 



M 



Here dR is the volume element of the three-dimensional 
rotation group given by 



dR = 4 sin 2 (-) du sin 6 dQ 

A simplified expression for BiPS in real space is 
(2L + 1 



(12) 



8tt 2 



dfl^ J dVtiy 2 C(hi,m) 
J dR X L (R)C(Rn 1 ,Rn 2 ) . (13) 



For statistical isotropic model condition = kqSlo can 
be recovered using orthonormality of x L (i?), 



X L (fl) X L '(fl)sin 2 -du; 
o 1 



?LL> 



The BiPS of CMB anisotropy computed from the maps 
measured by WMAP are consistent with SI, rulings out 
its radical violation [20J . An advantage of BiPS is that its 
rotational invariance allows for constraints to be placed 
on the presence of specific forms of CMB correlation pat- 
terns independent of the overall orientation in the sky. 



C. Reduced Bipolar coefficients- rBipoSH 

In order to extract information on the orientation of SI 
violation, or to detect correlation patterns in a specific 
direction in the sky, the Reduced Bipolar coefficients [23| , 
obtained as 



A 



LM 



oo L+li 

E E 

h=0l 3 =\L-h\ 



aLM 



(14) 



provide another set of measures. The summation of Bi- 
poSH over spherical wave-numbers l\ and l 2 , reduces the 
cosmic variance rendering these measurable from the sin- 
gle CMB sky map available. 

Note that the summation involves both the terms 
Afi, and Aff*, that are related via symmetry prop- 
erties Eq. Thus for any such combination where 
h + h — L is odd, these two terms will cancel each other 
leaving no contribution to the summation. The reduced 
Bipolar coefficients A^m, by definition have the following 
symmetry 

A LM = (-1) M A* L _ M , (15) 

which indicates Alq are always real. When SI condition 
is valid, the ensemble average of Alm vanishes for all 
non-zero values of L, 



< A LM >= 0, V L ± 0. 



(16) 



These Alm coefficients fluctuate about zero in any given 
CMB anisotropy map. Therefore, a statistically signifi- 
cant deviation from zero would confirm violation of SI. 
Unlike BiPS, reduced Bipolar coefficients are sensitive to 
orientation, hence they can assign directions to correla- 
tion patterns of the map. 



D. Bipolar map 

It is possible to visualize correlation patterns using the 
Bipolar map constructed from the reduced Bipolar coef- 
ficients Alm as [23[, 



6(n) 



LM 



(17) 



The Bipolar map from Alm is computed in the same 
way as the temperature anisotropy map from a given 
set of spherical harmonic coefficients, a/,„. Bipolar map 
can also be represented in terms of Tripolar Spherical 
Harmonics of zero angular momentum (see appendix [B] 
for details), 

9(n) = J2 I dn Al dn fl2 C(n 1 ,h 2 )(-l) h+h ^/(2L + l) 

L,h,l 2 

8xL{Y L {n) ® {Y h { ni ) ® Y h (n 2 )} A } 00 . (18) 
The tripolar spherical harmonics are expressed as (24| 
{Y h (hi) <E> {Y h (h 2 ) <S> Y h (n 3 )} l23 } LM = 



E ^mii23m23^'2m™/3m 3 ^l™l (ni)Yi 2 , n2 (h 2 )Yi 3m3 (n 3 ) , 



where the summation is carried over mi,m 2 ,m 3 , and 
m 23 . The transformations under rotations of tripolar 
spherical harmonics are identical to spherical harmonics. 
In particular, the tripolar scalar harmonics, which are 
invariant under rotations, can be expressed as follows, 

{Y h {hi) ® {Yi 2 (n 2 ) ® Y h (n 3 )}A}oo - (-l)' 1+ ' 2+ ' 3 <5 A h 

( h h h S 
ymi m 2 m 3j 



E (I 1 , ™~ ^ n ) Y hm 1 (hi)Yi 2m2 (h 2 )Yi 3m3 (h3). 



mim2m3 
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Orthogonality and normalization relation is as follows, 

dQ,f ll dSlf l2 dQ.f l3 {Y h (h 1 ) ® {Y h (h 2 ) ® li 3 (n 3 )}A}iM 
{^(ni) (8 {^(n 2 ) ® ^(n 3 )MW 

From Eq. (| its evident that under SI the Bipolar 
map is invariant under the rotations, since tripolar scalar 
harmonics are rotationally invariant and C(n\,h-z) = 
C(Rnx,Rn2). Hence, the map gets contribution only 
from the monopole term Aoo, 

Also, if the temperature map is rotated by a element 
of rotation group, "i?" then Bipolar map also rotates 
identically (see Appendix |B|) . For example, if you rotate 
the temperature map about the z-axis by some angle "a" , 

AT(R{9, 4>)) = ai m Y lm (0, <t> - a), 

lm 

the Bipolar map will also be rotated about z-axis through 
same angle "a" 

Q(R(6, <j>)) = A LM Y LM {0, <P - a). 

LM 

However, the Wigner-D matrices in the two cases will be 
different because of different m (or M) values. 

III. BIPOLAR REPRESENTATION OF CMB 
CORRELATION SYMMETRIES 

The homogeneity and isotropy of cosmic mi- 
crowave background points to the Friedmann-Robertson- 
Walker(FRW) model of universe. Flat FRW model ad- 
equately describes the observed local properties of the 
universe, but the fact that universe with same local ge- 
ometry can admit different global topology has been ap- 
preciated since the advent of post GR modern cosmology. 
This is because Einstein's equations describe local prop- 
erties of the spacetime and can only constrain, but not 
determine, the global topological structure. 

Symmetries of the space are preserved in the corre- 
lation function and global topology modifies correlation 
function. The simply connected (topologically trivial) 
hyperbolic 3-space R 3 , and the flat Euclidean 3-space £ 3 , 
are non-compact and have infinite volume. There are nu- 
merous theoretical motivations, however, to favor a spa- 
tially compact universe (j-u| ■ Compact topologies (more, 
generally, multiply connected space) break the statistical 
isotropy of CMB in characteristic patterns and induce 
a cutoff in the power spectrum because of finite spatial 
size 0, I25U27I ] . Theoretical possibilities include com- 
pact Euclidean and Hyperbolic 3-spaces which require 



the space to be multiply connected. The compact hyper- 
bolic manifolds are not globally homogeneous and they 
turn out to be not of much interest for the class of sym- 
metries considered under the scope of this paper. 

Simply connected universes are statistically isotropic, 
i.e. C(ni,n2) = Cin^.n^). In contrast, all compact 
universe models with Euclidean or hyperbolic geometry 
C '{hi ,712) are statistically anisotropic. The isotropy of 
space is broken in multi-connected models; this break- 
ing of symmetry may be apparent through the presence 
of some principal directions. In a cylinder, for instance, 
which is compact in one dimension and infinite in the 
other two, the metric tensor is exactly the same at ev- 
ery point hence it preserves local homogeneity. However, 
it is not globally isotropic and does not have the maxi- 
mal symmetry. It is noteworthy that globally anisotropic 
models do not contradict observations, since the homo- 
geneity of space and the local isotropy can ensure the 
observed isotropy of the CMB, however can influence 
the spectrum of density fluctuations. Multiply-connected 
models with zero or negative curvature can be compact 
in some, or all their dimensions. For instance a toroidal 
universe, despite its zero spatial curvature, has a finite 
volume which may in principle be measured. It contains 
a finite amount of matter. A cylindrical universe (in the 
sense that the spatial sections are cylinders), on the other 
hand, is noncompact in one dimension only and has an 
infinite volume, although with a finite circumference in 
the principal direction. 

Homogeneity and isotropy are experimentally con- 
firmed in the observations of distribution of luminous 
red galaxies 12811 . and the isotropy of CMB back- 
ground [1^, |30| |. Most of the studies in CMB as- 
sume statistical isotropy of the universe (FRW model). 
However, indications for a preferred direction in CMB, 
have motivated the study of departures from statistical 
isotropy These deviations can arise from non-trivial 
spatial topologies [1, [f| l8]-fl0ll, o r departures from the 
background FRW metric [IJ, l31j . Alternatively, statis- 
tical anisotropies might also arise from coherent mag- 
netic fields in the universe [1J, [l5|, [32j . Anisotropic Cos- 
mological models have been considered in the past and 
they lead to characteristic patterns in the CMB sky . 
The Bianchi template is an example of SI violation due 
to departure from background FRW geometry. Here we 
will discuss the signature of anisotropy due to existence 
of preferred axis (axes) on BipoSH. Such SI violations 
can arise due to non-trivial topologies as well as coherent 
magnetic fields. 

Since Bipolar formalism is sensitive to structures and 
patterns in the underlying two point correlation function, 
particularly the real space correlations, it is a novel tool 
to characterize statistical anisotropies [l7|, Ii9- 23] . Ro- 
tational symmetry about a preferred axis (say z) is the 
simplest way to break SI. 

In general, the correlation function may be decom- 
posed into isotropic and anisotropic parts [26|, 

C(ni,n 2 ) = C«(ni,n 2 ) + C^ni.na). (20) 
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where 

C^(nin 2 ) = C(n r n 2 ) = E C ; P ; (n r n 2 ) , (21) 

and the anisotropic part C^ A ' is orthogonal to the Leg- 
endre polynomials 

J dn hl J cM^ {fix ,n 2 )P l (h 1 -n 2 ) = . (22) 

This decomposition is useful in our study of the symme- 
tries of the CMB correlation patterns/structure that are 
explicit in real space. 

A. Statistical Isotropy (Rotational symmetry) 

Under SI, the correlation function is a function only 
of 9, the angle between the two directions, say, fi\ and 
fi2- Hence, C(ni,n 2 ) = C(fi\ ■ n 2 ) = C((9), and the cor- 
relation function can be expanded in terms of Legendre 
polynomials 

oc 

C (J '(fl) = -rY(2l + 1)CiPi(cos0), (23) 

47T ' 

1=2 

where Ci is the angular power spectrum. The summa- 
tion starts from 1 = 2, since I = and 1 = 1, respectively, 



monopole and dipole, are usually subtracted out. For SI 
the angular power spectrum C\ contains all the informa- 
tion. 

In Bipolar representation, the condition of SI for var- 
ious observables, described in the previous section, can 
be summarized as follows [13, [HI : 

. BipoSH : Affi = {-l) h C h (2h + iy/ 2 S hl J L0 6 M0 , 
• BiPS : k l = k 5 L q, 

. rBipoSH : A LM = E h (-l)' 1 ^ 1 (2^i + l) 1/2 '5LO^/o, 
. Bipolar map: 9 = I £ ; (-l)' y / IHp7J ; . 

Therefore, to test a CMB map for statistical isotropy, one 
should compute the BipoSH coefficients for the maps and 
look for non-zero BipoSH coefficients. Cosmic variance 
calculated for BipoSH under statistical isotropy is (see 
Appendix |C|) , 

4i{Kt)=C h C h [l + {-l) L 5 hh ] (24) 
for rBipoSH is, 

4i(Alm) = £ C h C l2 [l + (-1)^+^] (25) 
hh 

and for BiPS 0, 0, 



°sj(«l) 



l:2l>L 



E <t 2 2l + l +(-l) L (2£+l) + (l + 2(-l)V»]+]T E 4Q 2 1 Q 2 2 [(2L + 1)+F^ 2 ] 

" ' h h=\L-h\ 

+ 8 E^t"^[ E ^]» + i6(-i)^ E E ^ pa) 



i 2 =|i-ii| 

r 



l:2l>L 



where 

ii ( 3 L 

m\rri'2 — — l\ 77137714 — — Z3 M,M' — — L 

3 77I4Z1 771 1 — 7713/1 —777-2 



I3 —7713 ^17712/37774 



(27) 



and L is even. Statistically significant deviations from 
zero would mean violation of statistical isotropy. 



B. Cylindrical symmetry 

The correlation function must satisfy the symmetries 
of the underlying theory. In Friedman models the sym- 
metry group is SO (3), hence the correlation function is 



invariant under rotations; any breakdown of SI will re- 
duce this symmetry group. The simplest way to break SI 
is to introduce a favored direction in the sky, in such a 
case the reduced symmetry group is SO(2) or cylindrical 
symmetry. Assuming the favored axis to be z-axis, the 
rotational symmetry about z-axis for any arbitrary A(f> 
will require, 

cW(9i, <f>i, 2 , fa) = cM{p u <k. + N>, h, + A0). 

where n\ = (6±, <j}\) and n 2 = (# 2 , </> 2 ). The most general 
form of the correlation function in such a case is (see 
Appendix |D|) 

C (A) (0i,<£i, 02, <f>2) = E/m^^) co S m(^ - (Pa). (28) 

m 

Further, if the correlation function is invariant under the 
reflection, i.e., looking at a correlation pattern in the 
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sky one cannot distinguish whether we are looking up 
or down the preferred direction, then 

C(tt - X , fa, 7T - 2 , fa) = C{9 1 ,fa,9 2 ,fa), (29) 

which leads to a condition that h + h is even. BipoSH 
in such a case would be, or equivalently the covariance 
matrix will be [33J , 



- * c /~lhh 

<■ fl !imi a ! 2 m 2 > '— "mim2 lj | mi 



(30) 



where diagonal terms C l m of Cjj^i are called the cylindri- 
cal power spectrum, and \m\ > modes are the allowed 
frequencies for scale I. There may be correlations be- 
tween various scales called connectivity of fluctuations. 
The expression for Ci£* | in terms of the correlation func- 
tion is 



c hh = J_ (2h + l)(2l 2 + l)(h-m 1 y.(l 2 -m 1 y. 
|mi1 87r\/ (fi + mi)!(J 2 + mi)! 

P hmi (cos 6i)P hmi (cos 9 2 )f mi {9i,0 2 )d(cos 6>i)d(cos 9 2 ) 
Using Eq. (jDllj) in the appendix, we obtain 



AtX =[! + (-! 



V L ) °|m| U hml 2 -m - 



M0- 



When m is even, the functions P/ 71 will be odd or even 
functions of their arguments, depending on whether I is 
odd or even respectively. Similarly, for the odd m's. In 
both the cases when only one of l\ and l 2 is odd, the 
integral vanishes. Therefore we have to consider cases 
when both of them are either odd or even. In such a case 
l\ + 1 2 is even and hence Af 1 ^ vanishes for L = odd, 

Kh = KfJ^MO, where k = 0, 1, 2, 3, .... (31) 



.4 



LM 
hh 



Ai i SloSmoSi^ +Af J £5 M o 



LM ; 



(32) 



Therefore, the BipoSH present under cylindrical symme- 
try are A®® and Aff with even L. Using symmetry prop- 
erty of BipoSH {6$, under cylindrical symmetry we have 



LM 
hh 



AfJ^, i.e., the BipoSH are symmetric under the 



A 

exchange of and l 2 . There is another possibility that 
ai m 's have a gaussian distribution with different variance 
for each m mode corresponding to a particular I. This 
implies breakdown of SI, as power in each m mode is dif- 



ferent, C\t\ = (5/,/oC 



y |mi 

coefficients are, 



hh^\ mi \' 



and the corresponding Bipolar 



A hh = A^SLoSMoShh + Af^dMo^u 



1 LM 



LM i 



(33) 



In such a case non-zero BipoSH are A™ and Aff, where 
multipole moment is even and (L > 2). Furthermore, it 
may happen that a given model displays the degeneracy 



C'lmil = C/j and the rotational symmetry SO(3) of co- 
variance matrix is restored. The rBipoSH for cylindrical 
symmetry are, 

\ " aLM 



A 



LM 



hh 



AtX = AlmSl^Smo where k = 0, 1, 2, 3, .. 



Hence the Bipolar map for such a symmetry will be, 
9 (") = ^A LM Y LM {h) = ^A L0 8 Lak Y L0 {h) 

LM L 



(34) 



= J2\J^^ A ^L^kPL(cose) (35) 

Thus the map here looks like a sphere which is divided 
into latitude bands, or zones, without any longitudinal 
variation. 

A realistic example of cylindrical symmetry is a pri- 
mordial homogeneous magnetic field which breaks sta- 
tistical isotropy by inducing a preferred direction (e). 
Therefore, the correlation function between two points 
(n and n') depends not only on the angular separation 
between two points (n.n') but also on their orientation 
with respect to the magnetic field. This dependence of 
correlation function on angles between n and e (as well as 
n' and e) leads to correlation between I and / ± 2 modes 
[H| . The vector nature of the magnetic field induces off- 
diagonal correlations fl5) 

Di(m) =< a ; *_ lTO a/+i m >=< a,* +lm a/_i m > . (36) 

Here D; is the power spectrum of the off-diagonal ele- 
ments of the covariance matrix, and the correlation func- 
tion shows up as, 

<fl Iiim«! 2m2 > = S mim2 Si 1 i 2 Ci 1 + (37) 
u mim 2 

The BipoSH corresponding to this covariance matrix are 
0. 

Aff 2 = (-l) h (2h + l) 1/2 C h 5 h hSL SMo (38) 



AA^MO^C-ircS 



LO 
mh —m • 



The non-zero BipoSH in this case are A®® and Af^ i±2 . 

The reduced Bipolar spherical harmonic coefficients 
(rBipoSH) for this case are 

A LM = Y. S mo5lo{21 + 1) 1/2 C'i 
i 

+2 Y,(-±y n DiCt° lml+lm 6 M o. (39) 

in i 

These coefficients are non-zero only for li + 1 2 — L = even, 
thus L can take only even values. Finally, the Bipolar 
map is, 

6(n) = J2 Alo6l ^ Ylo ^^ k = 0,1,2,3 

L 

= -^Aoo + ^T A L0 S L . 2a P L (cos 9), a =1,2.. 
2 V 7r r 
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where A o and Alo are given by Eq. ([39)) . 

C. n-fold discrete Cylindrical symmetry 

Violation of SI also manifests itself in compact uni- 
verses with flat universal cover, which exhibits a n-fold 
rotational symmetry about an axis. There are six pos- 
sible compact models of the universe having a fiat uni- 
versal cover (UC) @ • These are visualized by identifying 
opposite sides of the fundamental polyhedra. The fun- 
damental polyhedron (FP) may be a parallelepiped. The 
possible identifications then are (figured]) 













F' 





1) E F G H - 
EFF'E' - 
EHH'F 

2JEFGH 
EFFE' 
EHH'E' 

3) E F G H 
E F F E' 
EHH'F 

4) E F G H 
EFFE' 
FGG'F 



E' F' G' H' 
HGG'H' 
FGG'F 

G' H' E' F 
HGG'H' 
FGG'F 

F G' H' E' 
HGG'H' 
FGG'F 

G' H' E' F 
- H'E'EM 
• G'H'HG 



FIG. 1: The locally Euclidean, closed, oriented 3-spaces. 

1) - opposite faces by translations. 

2) - opposite faces, one pair being rotated by angle tt. 

3) - opposite faces, one pair being rotated by tt/2. 

4) - opposite faces, all three pairs being rotated by n. 

The fundamental polyhedron can also be the interior of 
an hexagonal prism, with two possible identifications (fig- 
ure : 



1) E F G H I J <-»• G'H' I' J' E'F 
EE'J'J -«-»- GG'H'H 

I I' H' H -«-»- EFFF 

2) E F G H I J ■«-»- F G' H' I' J' E' 

EE'J'J ■*-»■ GG'H'H 
II' H' H -»-»- EFFF 



FIG. 2: The locally Euclidean, closed, oriented 3-spaces. 



1) - opposite faces, the top face being rotated by an 
angle 2ir/3 with respect to the bottom face. 



2) - opposite faces, the top face being rotated by an 
angle ir/3 with respect to the bottom face. 

Correlation function having a n-fold rotational symmetry 
about z-axis can be written as 

CW(9 U fa, 2) fa) = c^ A \e 1 ,fa + —,6 2 ,fa + —). 

n n 

This symmetry enforces mi + m 2 — nk, where n can 
be odd or even, depending upon the symmetry of the 
compact universe and k — 0,1,2,3... . Thus the general 
form of correlation function is (see Appendix IE)) . 

c (A) (0i, fa,e 2 , fa) = ]T f mitma (0i,h)x 

mi ,7712 

e^ m ^ +m ^S mi+m2 , nk , fc = 0,±l,±2... (41) 



Corresponding Bipolar spherical harmonic coefficients 
under the symmetry ea. (|40l) are, 



a lai _ i \- / (2/ 1 + i)(2i 2 + i)(/ 1 -m 1 )!(i 2 -T^)i" LM r 2 * r 2 * r r 

lh ~^hJ (h + mi m + m2 y. ^ 1+m2 , nk c hmihm2 j Q y o j o j Q 

C(9 1 ,fa,e 2 ,fa)e^ m ^ +m ^P hmi (cos6 1 )P l2 „^^ (42) 



All possible Euclidean models of compact universe ex- 
hibit reflection symmetry about the xy-plane. The cor- 
relation function under reflection symmetry is, 

C(6 1 ,fa,6 2 ,fa) = C(<K- 6 1 ,4> 1 ,k-6 2 , fa), (43) 

also under reflection of the coordinate system about x-y 
plane the spherical harmonics transform as, 

Y lm (n -6,fa = (-l) l+m Y lm (0, fa • (44) 



Therefore, reflection symmetry demands, 

P™ 1 (co S & 1 ) = P i 7 1 (co S (7r-^)) =P^(- Wa 1 ). (45) 

This implies l\ + mi is even and similarly l 2 + m 2 . Here 
we have used the symmetry property of Legendre poly- 
nomials, Pi m (— x) = (— l) l+m Pi, n (x). Interestingly, from 
symmetries of spherical harmonics, one can show that n- 
fold symmetries are ruled out for odd n (see Appendix 
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Topologically compact universes exhibits even fold 
symmetry, but the emergence of this fact from the sym- 
metry of two-point correlation pattern itself is neverthe- 
less instructive. Therefore, we need to look at the cases 
when n is even. 



D. Even- fold Cylindrical symmetry 

Even fold symmetry refers to the case when n is even. 
For compact topologies this is always the case, for in- 
stance, Dirichlet domain (DD) of a T 2 toroidal universe 
[35 | , and a T 3 have a 4-fold symmetry, that of a hexago- 
nal prism has a 6-fold symmetry and a squeezed torus has 



J 



2 fold symmetry. This symmetry puts another restriction 
on correlation function, 



C(0l,0l,0 2 ,&) =C(01,-0 



1,V2, 



0- 



(46) 



Hence most general correlation function under even-fold 
symmetry is (see Appendix [Ej) , 



mi+m2 ,nk 

mi ,??i2 

cos(mi<£i + m 2 ^ 2 ). 
Therefore Bipolar spherical harmonic coefficients are, 



^ = [i + (-i)' i+,a - i ] E 



'(2Zi + 1)(2J 2 + l)(Zi - mi )!(I a - m 2 )> iM 



(h +mi)!(/ 2 + m 2 )! 



< -^iimii 2 m2 0m i+ m 2,™fc 



27T />27T />7T />7T 



Jo Jo Jo 



C(0i , fa , 9 2 , fa ) cos(mi fa + m 2 </> 2 ) P™ 1 (cos 6>i ) P™ 2 (cos 2 ) dfa dfad{cos 9 X )d(cos 9 2 ) . (47) 



Reflection symmetry allows A; j on ly f° r even values of 
L. For odd indices coefficients vanish 



IV. BIPOLAR MAP: EXAMPLE OF THE 
BIANCHI TEMPLATE 



a LM aLMc X 
A hh — A hh°MnkOL2a, 



k = 0,1,2..., a = 0,1, 2, 3.. .(48) 



Here nk is even, therefore for an even- fold symmetry i.e., 
(n = 2, 4, 6, 8...) M is even and multipole moment L > 2 
picks up only even values. Example: for a 2 fold sym- 
metry possible values of M = 2,4,6,8...., similarly for 
4-fold symmetry M — 4,8,12,16... and so on. Another 
symmetry property of BipoSH is Eq. ©, so under even- 
fold rotational symmetry l\ + 1 2 + mi + m 2 = even which 
implies h+l 2 is even hence Af^ = Af^f i.e., BipoSH are 
symmetric under the exchange of l\ and l 2 . Note that for 
all possible compact flat spaces BipoSH vanish for odd 
indices, and the fact that two-point correlation function 
is invariant under reflection about xy-plane plays a piv- 
otal role in restricting non-zero BipoSH only to even Us. 
The rBiposh for even universe with even fold symmetry 
are, 



Alm — A LM S Mnk S L2a 



Bipolar map will be, 



= ^ A LM Y LM (n)5Mnk^L2a, 



(49) 



(50) 



LM 



L will take only even values and M will run from —L to 
L, subsequently picking up even values. 



Now we will consider a Bianchi template as an exam- 
ple to show how a Bipolar map looks like for a given 
temperature map. The choice of Friedmann-Robertson- 
Walker (FRW) model as a model of our universe was 
initially due to its simplicity, and later because of ob- 
servational evidence which strongly suggests universe to 
be homogeneous and isotropic at large scales. However, 
the presently observed isotropy may not necessarily hold 
in the past and the universe may have been anisotropic 
in its early stages and tends to FRW only later as it 
evolves. Bianchi models are the simplest examples which 
have the property to isotropise as they evolve in future. 
Bianchi classification contains 10 equivalent classes giv- 
ing generic description of a homogeneous and anisotropic 
cosmology [l2| . The most general Bianchi type which ad- 
mits FRW at late time are Vllh and IX. However, the 
type IX re-collapses after a finite time hence do not come 
arbitrarily close to isotropy. Spiral pattern are charac- 
teristic signatures of VIIq and VI In models [H, HE HH • 
Jaffe et. al. proposed Bianchi Vllh models as an ex- 
planation of WMAP anomalies. Since class Vllh models 
resembles a universe with vorticity and hence can lead to 
bounds on the universal rotation in cosmological (CMB) 
data [37l ]. They proposed correction for some anomalies 
in the first year maps from WMAP, however, introduc- 
ing such corrections induces other features like preferred 
direction and violation of SI. Pontzen et al. calculated 
various temperature and polarisation anisotropy patterns 
which may be formed in Bianchi cosmologies [38j . Ghosh 
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et. al. analyzed the temperature map for Bianchi Vllh where super-script B signifies Bianchi, and fi(9) and 
template [13(. Given the temperature map for Bianchi f 2 {&) are parameters of the model which should be cal- 
VI Ih template, here we see how a Bipolar map actually culated numerically [l^ . 
looks like. 

The temperature map for Bianchi Vllh template is of 
the form 



AT B (9, 0) = h (0) sin + f 2 (9) cos <£>, 



(51) 



BipoSH for Bianchi template are, 



.4 



LM 



{Wi 1 i 2 (0 1 ,9 2 )C^ lhl 5Ma + Xi 1 i 2 (9i 1 9 2 )C[ j 1 i ^ 112 _ 1 6m-2 



LM 



JO 



Y h iMMGtfh-^MQ + Z hl2 (9i,6 2 )CfXi 5 M2}d(cos6i)d(cos9 2 ). 



Therefore, rBiposh are, 
Alm 



where 



n7T 
{ W h i 2 (9!, 9 2 )C l L l A f ll2l 5 M Q + Xi 1 i 2 (9i,9 2 )C{ J l M ll2 _ 1 



F ili2 (^^ 2 )^ 2 _ 1 5 M o + ^ 1 i 2 (^^2)QXl^2}d(cOS^)d(cOS& 2 



LM 



/ (2Z 1 + l)(2Z 2 + l)(i 1 + l)!(Z 2 -l)! 
(4tt) 2 (Zi ^ 1)!(Z 2 + 1)! 



{Wi)h{e 2 ) + i{fi(9i)f 2 (0 2 ) - f 2 (9i)h(9 2 )) + 
f 2 (9 1 )f 2 (9 2 )}p- 1 (cos9 1 )P l \(cos9 2 ), 



(52) 



(53) 



/a (0i ) / 2 (0 2 ) jP,; 1 (cos ^ (COS 2 ) , 



*U = ^t^+i)^ - i)! + ^ ^^^ + + M'OAWO) + 

/2(0l)/2(02)}P ; 1 1 (cOSt9 1 )^7 1 (cOSc9 2 ), 



/a (^1 )/a (02)}^ (cos 0! )i£ (cos 2 ) . 



Hence, Bipolar map for a Bianchi template looks like (see ap- 

pendix [Gj , 

Alm = A LM S M ,k k = 0, ±2 . (54) 

Possible values of M are 0,±2. Note that Alm exists 
only for l\+l 2 — L = even, and vanishes otherwise. Keep- 
ing in mind the reality of two-point correlation function, = AloYlq^, <ft) + 2 Xl- 2 Gl{6) cos 2<ft 
i.e., A LM = (-1) M A L _ M , here we have A L2 = A* L _ 2 . l l 

Since Alm coefficients are complex numbers we can de- _ 2 Zl- 2 Gl{9) sin 2(f), (55) 

fine, X L m = $t(A LM ) and Zlm = %(A LM )- Therefore, L 
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FIG. 3: Temperature map for Bianchi VIIh- 



where 



GiW -(^\/ An(2L + l) (56) 

P L _ 2 (cos0) 2(21- + 1)P L (cos 0) P L+2 (cos6 



2L-1 (2L- l)(2L + 3) 2L + 3 







FIG. 4: Bipolar map for Bianchi VIIh- 



Thus, a spiral pattern in temperature map will show 
up as double spiral pattern in Bipolar map. 



V. CONCLUSION AND DISCUSSION 



azimuthal space. We found that the corresponding pat- 
tern in Bipolar space becomes a double spiral having a 
quadrupole in azimuthal space of the Bipolar map. 

Another application is in case of homogeneous isotropic 
models where an anisotropic topological identification 
has been imposed. As an example, if the space is com- 
pact in one (or more) direction(s), the statistical isotropy 
is broken due to introduction of preferred direction(s). 

We calculate BipoSH when this preferred direction is 
introduced. We have shown here that for compact topolo- 
gies, symmetry requirements can restrict BipoSH to even 
multipole moments, i.e., BipoSH vanish for odd indices 
for all kind of physically plausible models of flat multi- 
connected universe. Hyperbolic manifolds do not have 
the desired symmetry and hence we expect odd multi- 
poles to be non-zero in these manifolds. Hence, we have 
a tool to distinguish different topologies. In case of ho- 
mogeneous magnetic fields we have shown that BipoSH's 
are restricted to even L and M = 0. 

A new representation of Bipolar map has been pro- 
posed. Further work needs to be done in this direction 
to extract new information from this representation. 

This technique can be applied to polarization maps 
and it may prove to be a powerful method to decipher 
the topology of the universe, something on which general 
relativity is completely silent. The Bipolar formalism 
can also be applied to various anisotropic universes and 
can be used as a tool to distinguish various types of SI 
breakdown. 
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Representation of correlation function of CMB 
anisotropy in terms of Bipolar spherical harmonics pro- 
vides a novel approach to study violations of SI. Very re- 
cently the Bipolar representation has been used to quan- 
tify anomalies in the analysis of WMAP seven - year 
data [42j]. These anisotropics can arise due to departure 
from FRW metric (eg. Bianchi models), non-trivial spa- 
tial topologies (compact spaces) or from primordial mag- 
netic fields, among others. Here we have studied various 
measurable quantities of Bipolar formalism to quantify 
breakdown of SI. 

We studied anisotropic homogeneous cosmologies 
which leave a characteristic pattern on CMB. Like 
Bianchi VIIh temperature map which has a spiral pat- 
tern of a pair of cold and hot spots with a dipole in 



Appendix A: A review of topologically compact 
spaces 

Topologically compact spaces break the statistical 
isotropy, thereby introducing signatures in CMB corre- 
lation patterns. A compact cosmological model, Ai, is 
a Quotient space, constructed by identifying points of 
standard FRW space under the action of suitable dis- 
crete subgroup of motions r, of the full isometry group 
G of the FRW space. The isometry group G is the group 
of motions which preserves the distance between points. 
The simply connected infinite FRW spatial hypersurface 
with same constant curvature geometry is the universal 
cover (UC), Ai u , tiled by the copies of the compact space, 



11 



M. It can be spherical (S 3 ), Euclidean (£ 3 ) or hyper- 
bolic (H 3 )- The compact space for a given location of 
observer is represented as Dirichlet domain (DD), with 
the observer at its basepoint. Any point x of the com- 
pact space has an image, Xi — 7^x, in each copy of DD 
on the universal cover, where 7, e r. By construction 
DD represents the compact space as a convex polyhedron 
with even number of faces identified pair wise under _T. In 
cosmology, DD around the observer represents the uni- 
verse as seen by the observer and the symmetries of the 
correlation function are nothing but the symmetries of 
the corresponding DD. 

Correlation function of a scalar field, <P, on a compact 
manifold, M, can be expressed as (4pj |. 

rrii 

££(x,x') =^^P # (fc i )%(x)^*.(xO, (Al) 

i j = l 

where 

(V 2 + fc 2 )*y = 0. (A2) 

are orthonormal set of eigenfunctions of Laplacian 
on the hypersurface, having positive and discrete set of 
eigenvalues {kf} (fc 2 = and kf < kf +1 ) with multiplic- 
ities rrii. On a compact manifold, M, the set of eigen- 
functions and eigenvalues are not always easy to obtain 
in closed form (even numerically, for compact hyperbolic 
spaces). On the other hand, eigenfunctions ^"(/^x) of 
the universal cover (UC), A4 U , of a compact manifold 
usually known because of their simplicity (e.g., W 3 , S 3 
and £ 3 ), hence they can be used to compute the corre- 
lation functions £$(x, x') on the UC. For flat and hy- 
perbolic UC's the set of eigenvalues are continuous. The 
function P${ki) is the rms amplitude of the eigenmode 
expansion of the field <P, whose information lies in the 
physical mechanism responsible for the generation of <S>. 
The regularized method of images [25j , describes how cor- 
relation function on the compact manifold can be calcu- 
lated once the correlation function on the universal cover 
is known [lj| Hg], which is expressed as, 

££(x,x') = £a(x,x'). (A3) 

-yer 

This implies that correlation function on a compact 
space, M, can be expressed as sum over the correlation 
function on its universal covering space, M u , calculated 
between x and the images 7x^7 e r) of x'. The local 
homogeneity and isotropy demands that the correlation 
function on the UC is only a function of the distance 
between two points x and x' i.e r = d(x, x'). The corre- 
lation function on a compact universe with flat UC is, 

^(x,xO=Wfp^)^, (A4) 

Here, P&(k) can be determined from the early universe 
physical mechanism and di is the distance between the 



images of x and x'(do is the distance between original 
points). Summation implies summing over all images. 
Hence, correlation function depends not only on the dis- 
tance between two points and the distance of their im- 
ages but symmetry defines both the pair to have identical 
distance from their images i.e., take any two points on 
the last scattering surface and their corresponding images 
about xy-plane, correlation function will turn out to be 
invariant under this reflection. Figure (|5|) illustrates this 




FIG. 5: Images of two-pairs of point which are mirror images 
of each other about the XY plane in a T 3 space ( one dimension 
suppressed). 

point for a T 3 universe. The DD of a squeezed torus is 
shown in figure ^ . The choice of the axes here is a little 
bit more non-trivial. The xy-plane is not parallel to any 
of the faces of the DD or the FP, but still it would cut the 
LSS into two halves in such a way that there will be sym- 
metry under reflection about the xy-plane and on the xy- 
plane there will be 2-fold rotational symmetry. However, 




b = 90 - a 



FIG. 6: Dirichlet domain of a squeezed torus. 
here we point out that reflection symmetry does not hold 
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good for the compact spaces for which the opposite faces 
are glued together with a twist |39f . 

Topology of the universe leaves characteristic signa- 
tures on CMB. If the universe is finite and smaller than 
the distance to the last scattering surface (LSS), then the 
signature of the topology of the universe is imprinted on 
the CMB. For such a small universe LSS can wrap around 
the universe and will self-intersect. The intersection of 
the LSS, which is a 2-sphere with itself is a circle that 
will appear twice in the cosmic microwave background. 
Hence, there might exist pairs of circles which share cor- 
related patterns of temperature fluctuations. This circles 
in the sky [4l| method is a powerful and direct probe 
for detecting non-trivial spatial topology. The correlated 
patterns would be matching perfectly if the temperature 
fluctuation did not depend on the direction of observation 
and if the patterns were not distorted. However, the ob- 
served temperature fluctuations has direction dependent 
components, i.e. the Doppler effect and the integrated 
Sachs- Wolfe effect. Also observationally, galaxy cut and 
foreground removals can also distort the matching. How- 
ever, one can search for such patterns in CMB correlation 
function statistically. In a multi-connected space, there 
exist preferred direction(s) so that global isotropy is bro- 
ken. The angular correlation will then depend on two 
directions of observations and can also depend on the po- 
sition of the observer. This induces correlations between 
a; m 's of different I and m. Thus, another indirect probe 
is to search such patterns or signatures in the statistics 
of CMB temperature fluctuations. 



where ( ft ' 2 | are Wigner-3i symbols and are re- 
\mi M m 2 J 

lated to Clebsch Gordan coefficients in the following way, 

h imihm2 - VM3 + 1 \m 1 m 2 -m 3 ) 

Hence Bipolar map can be represented in terms of tripo- 
lar scalar spherical harmonics, 

0(h) = /"^ 1 dfla 2 C(n 1 ,n 2 )(-l)' 1+ ' 2 \/2TTl 

L,h,h 

Sx L {Y h (h x ) ® {Y L (h) ® Y h (h 2 )}x}oo (B4) 

The representation of the Bipolar map in terms of tripo- 
lar harmonic function makes the transformation proper- 
ties of the bipolar map under rotations explicit. 
In a rotated sky map, 

©'(-)= E /[ E 4fT*«-(fi) 

LMhh V X ,V^V,M' M" 

{^(ni)®Y^(n 2 )}i'M''] x 
{Y h (hi) ® Y* 2 (n 2 )}£ M Y iM (n)dQ fil dfi ft2 (B5) 

Using orthogonality of Bipolar spherical harmonics 

©'(«)= E E A vfT, D M"M,(R) 

{> M" 

YlmSi^Si^Sll'Smm' (B6) 



Appendix B: Bipolar map representation in terms of 
tripolar spherical harmonics 



Bipolar map is defined as 

e(h) = J2 A LMY LM (r 



(Bl) 



LM 
LM 



where A LM = J2hi 2 A hh> therefore 



LM hl 2 

Now using the expansion of Af^ we get 

0(h) = EE / d ^ Al / dft^C^!,^) x 

LM hl 2 

{Y h (hi) ® Y h (h 2 )}* LM Y LM (h) 
tripolar scalar spherical harmonics are defined as 

{Yi^hi) ® {Y L (h) ® Y h (h 2 )} x } m = (B3) 

1 M m 2 ) Ylimi ^ YLM ^ Yhm2 ^ 



^ E ( TO \ 



mi HI2M 



0'(h) = Y, A LM >Y D MM>(R)YLM(h) = 0(Rh)(B7) 



LM' M 



Thus when correlation pattern is rotated by "i?" , Bipolar 
map also rotates by "iZ" . 



Appendix C: Cosmic Variance of Bipolar Quantities 

Cosmic variance is defined as the variance of estimator 
of an observable constructed from a single sky map. In 
particular for BipoSH 



a*(Af£) =< (Atxr > < KZ > • (ci) 

Using Gaussianity of AT, one can analytically compute 
the variance of A^¥ . 



TLM _ * (_-,\m 2 r<LM 



I2— m 2 



(C2) 



therefore 



XLM 7*LM 



E E < a hnii a *l2m2 a hm' 1 a l2m' 2 > 



( 1 \m 2 +m' 2 t~iLM r~iLM 1 ' 
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Considering temperature field to be a Gaussian random 
field, one can expand the four-point correlation function 
in terms of two-point correlation function. Further, con- 
sidering the fact that under statistical isotropy the covari- 
ance matrix is diagonal eq.([3]) the above equation reduces 
to 

< Af^AhuJ > = C h C h 6 llh (2h + 1)S L0 6 M0 

+ C h C l2 [l + (-l) L S hl2 ] (C4) 

Also, we have 

< Af™ >= (2h + l) 1/2 C h 6 hl2 6 L() 6 M o (C5) 
Hence the cosmic variance is 

vlAKT.) = GiMl + (-l) L 6 hl2 ] (C6) 
Similarly for rBipoSH, 

\h+l 2 -L] 



therefore 



(C7) 



2kn 



m 1 +m 2 = —, fe = 0,±l,±2.... (D4) 

lAtp 

for zero fold symmetry ra\ + m 2 = which means mi — 
— m 2 , hence 

C^(htn 2 ) = J2 ^C^ lhm2 Y hmi (ni) 

li ,l 2 ,L,M,mi,m 2 

using the expansion of spherical harmonics in terms of 
associated Legendre polynomials, 



Y, 



I 111 



it) ^ - P ^ / (^ + i)G-w)!" Dm/ „„„^ 



therefore correlation function will be, 



(D6) 



C^Hh^) = £/ ro (0i,0 2 )e iro( ^ 2) (D7) 



Appendix D: Correlation function for Cylindrical 
symmetry 

Expansion of correlation function in terms of Bipolar 
spherical harmonics is, 

CW(n 1 ,n 2 )= £ Afjf£«^ 

Y hmi (hi) Yi 2m2 {fi2) (Dl) 

Now rotational symmetry about z-axis for any arbitrary 
A<j) implies, 

C {A) (0 1 ,(f> 1 ,6 2 , <fo) = C\6x,<h. + A0, 2 , cf> 2 + A0XD2) 
Therefore 

^ AfJ%Ci$% lla7na Y limi (9l, fa) Yi 2ni2 (0 2 , fa) 

1 1 , 12 , L , M , m 1 ,m.2 

E A if^ ll2m2 Y hmi (ei,fa + Afa 
Y l2m2 (d 2 ,fa+Acf>) 



where 



which means 



fm(6l,0 2 ) — — £ A hh C hrnl 2 -rn S MO 



h,h,L 



(2h + l)(2/ 2 + l)(Zi - m)!(Z 2 + m)\ 
V ii + m)!(i 2 - m)! 

P/ 1 "(cos0i)P i ; m (cos0 2 ) (D8) 

Symmetry ensures, 

C (A) (^.0i^2 ! 02) = C A (0 1 ,-0 1 ,0 2! -0 2 ) (D9) 
Imposing this symmetry we get, 

C (A) (0i, 0i, 2 , 0a) = £ fm(0i,02) cos m(0i - 02CP1O) 



This is the most general correlation function under zero 
fold rotational symmetry. Using 

Aff 2 = I <ff2 ftl dnft s C(n 1 ,n a ){ll 1 (fli)®y, a (fi2)}i M 



«(mi+mj)A^ 



(D3) and eq. ([DT0|) we get, 



a lm _ h , / , Vi+i 2 -Li 1 w V ( 2 ^i + !)(^2 + l)(/i - mggg - my. ^ LM 



2 -m*M0 



■/() 



(47r) 2 (fi+m)!(J 2 +m) 

P ! m (COS 01 )P 2m (COS 2 ) f m (&1 , 2 )d(cOS 01 ) d(cOS 2 ) . 



(Dll) 
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Appendix E: n-fold cylindrical symmetry 



Correlation function in such a case is, 

2~- 27T 

c (A > (e 1 ,<h.,e»,ih) = c {A) (0i , 0i + — , 2 , <t> 2 + — ) 

n n 



Hence, e^ mi+TO2 ^ » = 1, which implies mi + m2 
nfc,fc = 0,±1,±2,±3 

Most general form of correlation function will be, 



(El) 



This implies, 

^ A hh C hm 1 l 2 m 2 Y hm 1 (Sl, 0l)^ 2 m 2 (02, 02 ) = 

y2 A^hC^ i Y hmi (9 1 ,(f> 1 H )Y hm2 (e 2 ,(f>2 H ) 

I 

where 



C (A) (0i,0i,0 2 ,0 2 ) = 

J miTti2 v^l ' 2/ %i+m2,nfc 



f rfl fl^i - 1 >T / (gZi + 1)(2Z 2 + l)(/i - mi)!(f 2 ~ m 2 )! £M LM pmi,.„ fl \ P ™2 /___/, n 
/ mi m 2 (0i,0 2 ) - t- 2^ \/ ^ r^WfTX^Tw A hh C hmihm2 F l 1 (cosff 1 )F h (COS0 2 ) 



hl 2 LM 



[h + mi)\(l 2 + m 2 )\ 
I 



Demanding explicitly the two fold symmetry that holds 
for all even-fold symmetry, 

C(d u 0i, 2 , 2 ) = C(6 X , -0i, 2 , -0 2 ) (E2) 

this symmetry rules out the presence of sine terms in 
correlation function. Hence for even-fold symmetry cor- 
relation function reduces to, 

C ( ^(0i,0i,02,0 2 ) = fmum 2 (0l,02) 

S mi +m 2 ,nk cos(mi0i + m 2 2 ) (E3) 



Appendix F: Absence of odd- fold symmetries for 
compact spaces 

Let us take a compact space with reflection symmetry. 
It would demand that 

£7(tT - 01, 01,7T- 02, 02) =£7(01, 01,02,02). (Fl) 

It can be shown that 

£7(0!, 0i,02,0 2 ) = J2 (- i r 2 < «imi°r a m, > >< 

C^ lh - m A Y h(ni) ® il 2 (n 2 )} iM .(F2) 
The symmetry of spherical harmonics would ensure that 

{Yfe (7T - X , 0i) ® Y^TT - 02, 2 )}iAf 

= (-l) mi+m2 {l; a (01, 0i) "8 Y 2 {6 2 , 2 )} LM • (F3) 

This put together with equations (|F1[) and (|F2I) indicates 
that mi must be an even number. Since they are dummy 
indices, m-2 would be even too. Now let us consider such 
a space with n fold symmetry. Evidently 

C7(0i, 01 + — ,02,02 + — ) =C(01, 01, 02,02) (F4) 

n n 



since we know 

{1^(01,01 + — ) ® H 2 (0 2 ,02 + —)}LM = 

n ~ n 

27r(mi+m 2 ) ^ ^ _ 

n 

For Eq. (|F4j) to hold (mi + m 2 )/n must be even. Since 
mi and m.2 are even, n has to be even too. 



Appendix G: Bianchi template 

The temperature map for Bianchi template is written 

as 

AT(0,0) = /i(0)sin0 + / 2 (0)cos0. (Gl) 
Bipolar map can be expressed as, 

Q ( fl ) = EE E f dn^dfl^ < AT(ni)AT(n 2 ) > 

LM l±l 2 m 1 m 2 

C , "n 1 i 2 m 2 ^ mi («l)^m 2 («2)>LAf(n). 

The integrals over contribute only for m = ±1 other- 
wise it vanishes. Constraint on the values of mi = ±1 
and m-2 = ±1, admits only M = 0, ±2. Reduced Bipolar 
coefficient is then 

A lm =J2Y, fd^cM^ < AT(h 1 )AT(h 2 ) > 

hl 2 mim 2 

Alm exists only for l\ + l 2 — L = even, and vanishes 
otherwise and the reality condition demands Alm, he., 
Alm = (— l) M ^-ji,-M- Now Bipolar map is 



e{n)=Y,A LM Y LM {h) 



(G2) 



LM 
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but for Bianchi template it will be 

6(n) = J2 ^LoYLo(n) + Y A L iY Li (h) 



L>2 



L>2 



which can be written as 

6(n) - Y, A LoY LO (n) + Y A *L-2Y£_ 2 (fi) 



L>2 



Y A L-2Y L ^ 2 (n) (G3) 



L>2 

Since A' LM s are complex numbers, we define 
Alm = Xlm + iZ L M and A* LM = X LM - iZ LM (G4) 



and the Bipolar map (|G3[) can then be written as 

9(h) =Y A loYlo(ti) + Y X L-2( Y l-2(n)+Y L _ 2 (n)) 



+ iYZL-2( Y L-2(n)^Yl_ 2 (h)). 

L 



Defining, 
G L (9) = 



+!)(£ + 2) 



(sin6») 2 y 4tt(2L+1) 
P l _ 2 (cos0) 2(2L+l)P L (cos#) P L+2 (cos0) 



2L-1 (2L- l)(2L + 3) 2L + 3 

the Bipolar map is represented as 

Q(6, (b)=Y A loYlo(0, cj>)+Y X L -2G L (9)2 cos 20 

L L 

-^Z L _ 2 G L (0)2sin2$G5) 

L 

where we have used expansion of spherical harmonics in 
terms of associated Legendre polynomials 



Yt 



I 111 



and 



J±24> 



Y, 



l±2{ 



(L - l)L(L + 1)(L + 2) 



(sin6») 2 y 4tt(2£ + 1) 

P L - 2 (cos6>) 2(2i + l)P L (cos6») P L+2 (cos( 



. 2L-1 (2L- l)(2L + 3) 2L + 3 

Appendix H: Useful Mathematical Relations 

Orthonormality of spherical harmonics 



(2ftftF Jimi (ft)Y;* OT2 (n) = ^i 2 i5 miro2 (HI) 
Symmetry property of spherical harmonics 

Y;jh) = (-l) m Y^ m (h). (H2) 



Spherical harmonic expansion of Legendre polynomials 



A-7T 

Pi(n ■ n') = E ijm(n)^m(ft')- (H3) 



Property of legendre polynomial 



(H4) 



(Z + m)!' 1 

Symmetry properties of Clebsch-Gordan coefficients 
C^/j = (-l)° +6 - C C h 7 aa , (H5) 

/^fC7 / i \a-\-b— C(~iC— 7 

Summation rules of Clebsch-Gordan coefficients 
Y C Zbf3 C aab/3 = S cc >5 7Y {abc}{abc'} 

a/3 

Y. C alb(s c Zb>i3> = ^^j 6 bb'S^{abc}{ab'c} 

ay 

Y C alb0 C Z>b0> = $ aa >S{H3>{abc} 

C~f 



where 



Y[ = [(2a+l)(26 + l)....(2c+l)] 1 



/2 



(H7) 



and where {abc} is defined by 

r l \ _ / 1 if a + 6 + c is an integer , , 

^ •' 1 otherwise ^ ' 

and where a, b and c satisfy triangle inequality \a — b\ < 
c < (a + b). 

Tripolar spherical harmonics are expressed as, 
{Y h {hx) ® {Y h (h 2 ) <E> Y h (h 3 )} l23 } LM = 

/] Gl 1 m 1 l23m2 3 Gi2m 2 l^m 3 Yl inii (^l)F 2 m 2 (^2)F 3 m 3 ( 
mim2in 3 m23 

Tripolar scalar spherical harmonics are defined as 

{^(ni) ® {F L (ri) ® F 2 (n 2 )} A }oo = (-l)' 1+ ' 2+i (H9) 
li L Z 2 



^ E (mi 



M m 2 



^imi(ni)Yw(n)Yj 2m2 (n 2 ) • 



where f 1 5l ^ 2 | are Wigner-3i symbols. 
V mi M m 2 y 

Orthogonality of tripolar spherical harmonics, is given as 

dn^d^tffi^YjJfti) ® {Yi 2 (n 2 ) (8 Y h (h 3 )}x} LA 

{^(ni)®{^(n 2 )®^(n 3 )h'}W' 
= Si^Si^Si^Sw'Sll'Smm 1 ■ 
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